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ABSTRACT 


The  nerve  agent  VX  (O-ethyl  S-[2-(diisopropylamino)ethyl]methylphosphonothiolate)  has  been 
shown  to  decompose  in  contact  with  concrete  surfaces.  We  have  previously  shown  that  31P  NMR 
can  be  employed  to  study  1)  the  adsoiption  of  VX  into  concrete;  2)  the  decomposition  kinetics;  and 
3)  detect  decomposition  products.  Recently  reported  results  show  a  dramatic  effect  of  the  VX 
decomposition  kinetics  on  droplet  size,  with  0.01  pL  droplets  reacting  within  a  few  hours  in 
contrast  to  months  observed  for  VX  drop  sizes  on  the  order  of  several  pL.  The  objective  of  this 
paper  is  to  probe  the  effects  of  temperature,  moisture,  concrete  source  and  droplet  size  on  the  rate  of 
hydrolysis.  In  the  present  study,  31P  NMR  showed  that  the  hydrolysis  of  nerve  agent  simulant 
DMMP  (dimethyl  methylphosphonate)  on  concrete  increased  with  increasing  temperature,  amount 
of  water  added,  and  pH  of  the  concrete  used.  No  effect  on  hydrolysis  rate  as  a  function  of  droplet 
size  in  the  0.05  to  2.0  pL  range  was  detected. 


INTRODUCTION 

The  rate  of  decomposition  of  chemical  warfare  agents  on  substrates  commonly  present  in  a 
battlefield  or  urban  environment,  such  as  concrete  and  sand,  is  important  since  the  decomposition  to 
non-toxic  products  will  determine  when  the  area  will  be  safe  to  reenter  without  protective  gear. 

The  nerve  agent  VX  (O-ethyl  S-[2-(diisopropylamino)ethyl]methylphosphonothiolate)  has  been 
shown  to  decompose  in  contact  with  concrete  surfaces.1  This  decomposition  was  shown  to  proceed 
with  an  initial  hydrolysis  of  12%  of  the  VX  that  had  a  half-life  of  2.2  hours,  followed  by  a 
subsequent  decomposition  that  had  a  half-life  of  28  days  to  3  months  on  crushed  and  monolithic 
concrete,  respectively.2  Recently  reported  results  show  a  dramatic  effect  of  the  VX  decomposition 
kinetics  on  droplet  size,  with  0.01  pL  droplets  reacting  within  a  few  hours3  in  contrast  to  months 
observed  for  VX  drop  sizes  on  the  order  of  several  pL.  Differences  between  these  experiments 
were  the  source  and  pH  of  the  concrete,  drop  size,  neat  versus  dissolved  VX,  and  crushed  versus 
chunk  concrete.  This  work  involves  a  systematic  study  of  the  variables  of  time,  temperature,  pH, 
moisture,  concrete  source,  and  drop  size  on  the  hydrolysis  of  the  simulant  DMMP.  The  data 
obtained  from  these  studies  will  be  used  as  a  guide  for  future  experiments  with  VX. 

EXPERIMENTAL 

MATERIALS.  Dimethyl  methylphosphonate,  (DMMP),  CAS#  Registry  Number  756-79-6  was 
purchased  from  Aldrich  Chemical  Company  and  used  as  received.  31P  NMR  of  the  starting 
materials  indicated  that  it  was  the  correct  compound.  Concrete  samples  8GWW,  CB1-27C,  and 
CB1-37C  were  all  made  using  sieved  Quikrete®  concrete  that  had  been  sieved  through  a  10-mesh 
screen  to  remove  all  of  the  large  aggregate.  The  concrete:  water  ratio  (by  weight)  was  4.3:1  for 
concrete  CB1-27C  and  4.5:1  for  concrete  CB1-37C.  One  method  of  characterizing  the  samples  was 
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to  put  the  concrete  in  water  for  1-3  hours  and  measure  the  pH;  8GWW  and  CB1-37C  gave  a  pH  of 
8-9,  whereas  CB1-27C  gave  a  pH  of  9.5  to  1 1.  Concrete  CZ  was  received  from  the  Czech  Republic 
(the  aggregate  in  these  samples  was  less  than  2mm  diameter,  and  its  pH  was  ~1 1). 

METHODS.  The  concrete  chunks  were  inoculated  with  DMMP  and  then  with  water,  which  was 
put  directly  on  top  of  the  DMMP.  In  general,  a  total  of  4  pL  DMMP  were  used,  with  8,  24,  or  40 
pL  water,  which  gave  molar  ratios  of  12,  36  or  60  respectively.  Thus,  all  of  the  kinetic 
measurements  were  pseudo-first-order  in  [H20].  The  sample  was  put  into  a  1 0mm  NMR  tube,  and 
the  initial 3 1  P  { 1 1 1  j  NMR  total  integrated  peak  area  was  measured  using  a  1 0mm  liquids  probe  on  a 
Varian  Inova  wide  bore  instrument.  After  allowing  an  estimated  half-life  of  time  to  pass,  the 
3 1 P  J 1 1 1}  NMR  total  integrated  peak  area  was  remeasured;  it  was  not  possible  to  obtain  resolution 
between  reactant  and  product,  since  the  peaks  were  about  7,000  Hz  wide.  The  overall  signal 
intensities  indicated  that  at  least  80%  of  the  original  31P  nuclei  were  detected.  The  sample  was 
then  covered  with  ca.  1  mL  water,  and  allowed  to  sit  for  30  to  180  minutes  to  allow  the  compounds 
to  diffuse  from  the  concrete.  The  solution  peaks  were  20  to  50  Hz  wide,  and  thus  the  individual 
components  could  be  identified,  using  ^Pl'H},  'H  and  31P  NMR.  31P{'H}  NMR  total  integrated 
peak  area  was  measured  for  all  components  and  compared  to  the  original;  all  of  the  31P  present 
initially  was  found  in  the  aqueous  solution.  Rate  constants  were  derived  from  the  %DMMP 
remaining  (the  reaction  in  solution  was  shown  to  be  first-order  in  [DMMP];  it  was  assumed  that  the 
reaction  on  the  concrete  is  also  first  order.  Only  reactions  that  had  1 0%  to  90%  DMMP  present 
were  used  for  these  calculations.  Chemical  shifts  were  measured  relative  to  an  external  phosphoric 
acid  standard. 

Statistics.  JMP  Statistical  Discovery  Software,  Professional  Edition,  by  SAS  Institute  Inc.  was  used 
for  the  experimental  design  and  data  analysis. 


RESULTS  AND  DISCUSSION 

STABILITY  OF  DMMP  ON  AN  AMBIENT  CONCRETE  CHUNK 

The  31P  NMR  spectrum  of  the  DMMP  on  ambient  concrete  was  measured  as  a  function  of 
time;  the  parameters  observed  were  the  31P  relaxation  time,  Ti,  total  integrated  peak  area,  line  width 
and  chemical  shift. 

The  Ti  was  largely  a  function  of  the  concrete  used,  and  no  trend  with  time  was  seen  (Figure 
la).  The  Ti  was  less  than  0.1  seconds  for  concrete  G8WW,  0.5  to  2  seconds  for  concrete  CB1-27C, 
and  0.1  to  0.5  second  for  concrete  CB1-37C.  Possible  reasons  for  the  short  Ti  in  concrete  versus 
the  12s  Ti  in  solution  are  paramagnetic  impurities  and/or  restricted  mobilities  in  the  concrete.  Ti 
was  also  measured  as  a  function  of  drop  size  (Figure  lb).  Statistical  analysis  of  the  groupings  of 
drop  size  using  the  Student  t-test  at  a=0.05  showed  that  the  0. 1  and  0.25  pL  samples  were 
significantly  different  from  the  0.5  and  2  pL  samples. 

The  total  integrated  peak  area  of  the  DMMP  shifted  slightly  based  on  the  exact  amount  of 
DMMP  that  was  added  to  the  sample.  Thus,  the  average  of  all  of  the  integrated  peak  areas  for  that 
sample  was  taken,  and  each  individual  integrated  peak  area  was  normalized  by  the  average,  to  give 
a  normalized  intensity.  If  there  were  no  effect  of  sample  age  on  the  total  integrated  peak  area,  the 
normalized  integrated  peak  areas  would  be  randomly  distributed  around  the  1.0  horizontal  line 
(Figure  2a).  The  expected  range  would  be  -20%,  which  represents  the  repeatability  of  this 
experiment.  The  three  samples  that  were  examined  in  detail  all  have  intensities  within  the 
anticipated  20%  range,  but  the  intensity  gradually  diminishes  with  time.  Possible  explanations  are 
a  slow  evaporation  of  the  DMMP  or  a  product  that  has  an  extremely  broad  line  width  and/or  a  very 
longTi. 

The  line  widths  of  DMMP  on  a  number  of  concrete  samples  run  in  the  liquid  probe  were 
measured.  The  range  was  4,000  to  10,000  Hz;  most  peaks  were  6,500  to  8,500  Hz  wide  (Figure 


2b).  This  broadness  precluded  the  resolution  of  DMMP  and  its  hydrolysis  product,  MMP  (Figure 

3) .  The  chemical  shifts  observed  ranged  from  34  to  46  ppm,  regardless  of  the  concrete  used  (Figure 

4) .  By  comparison,  the  chemical  shift  of  DMMP  in  distilled  water  and  a  variety  of  aqueous  buffers 
ranged  from  38  to  41  ppm. 


Droplet  Size 


Figure  1.  a)  Plot  of  relaxation  time,  Tl5  versus  the  sample  age.  Code:  D=93G,  ■=CB1-27C, 
♦=8GWW,  A=CB1-37C.  b)  Plot  of  Tj  versus  Droplet  Size  for  4pL  DMMP  on  concrete  CB1- 
27C. 
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Figure  2.  a)  Normalized  Intensity  versus  time  for  select  samples  of  DMMP  on  concrete  CB1-27C: 
□=93G,  +=13 1A,  4<=77B.  b)  Plot  of  the  peak  line  width  versus  sample  age.  Code:  IH=93G, 
■=CB1-27C,  ♦=8GWW,  A=CB1-37C. 
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Figure  3.  Chemical  shift  of  DMMP  on  Various  concrete  samples  versus  sample  age  .  Code: 
□=93G,  ■=CB1-27C,  ♦=8GWW,  A=CB1-37C. 
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Figure  4.  Chemical  shift  evolution  of  DMMP  on  concrete  chunk  CB1-27C  with  added  water. 


KINETICS  OF  DMMP  HYDROLYSIS  ON  A  MOIST  CONCRETE  CHUNK 

The  rate  of  DMMP  hydrolysis  as  a  function  of  added  water  and  temperature  was 
measured  by  incubating  the  DMMP  and  water  on  the  concrete,  then  removing  the  reactant  and 
product  from  the  concrete  with  water.  The  extent  of  the  reaction  was  calculated  from  the  relative 
%DMMP  remaining,  as  measured  by  31P  NMR.  The  water  and  temperature  parameters  were 
varied  according  to  an  experimental  design.  The  data  were  fit  to  the  equation: 

Log  ti/2  =  a  +  b*uL  water  +  c*temperature. 

The  model  was  accepted  if  it  had  R2  >  0.9;  significant  factors  were  chosen  as  those  with  P<0.05. 

An  acceptable  model  (Figure  5)  was  found  for  concrete  CB1-37C  (R2=0.92): 

Log10  ti/2  =  4.0  -  (0.014*pL  water)  -  (0.038*temperature). 

The  standard  errors  of  the  coefficients  were  0.2,  0.003  and  0.004  respectively. 

The  data  were  analyzed  using  the  logio  of  the  half-life  because  the  range  of  the  data  is  large  and 
because  the  time  periods  of  interest  are  minutes,  hours  and  days  rather  than  one  hour  versus  ten 
hours.  Thus,  the  data  are  sufficient  for  order  of  magnitude  predictions,  but  cannot  precisely  predict 
if  a  reaction  will  take  20  or  40  hours.  When  the  calculated  and  observed  logio(ti/2)  were  compared, 
the  ratio  ranged  from  0.5  to  2.0,  with  a  mean  of  1.0.  Some  of  the  variability  in  the  replicate 
measurements,  (which  was  also  reflected  in  the  standard  errors  of  the  coefficients  in  the  equation) 
was  due  to  the  fact  that  concrete  is  a  heterogeneous  material.  Although  much  of  the  aggregate  was 
removed  from  concrete  CB1-37C  before  use,  since  the  samples  were  small,  heterogeneity  was  a 
likely  factor  in  the  variability  of  the  half-lives. 
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Figure  5.  Plot  of  temperature  vs.  pL  water  ,  the  contours  are  the  half-lives,  in  hours,  for  DMMP 
hydrolysis  on  concrete  CB1-37C.  The  rectangles  show  the  conditions  under  which  the  data  was 
collected. 


EFFECT  OF  DMMP  DROPLET  SIZE  ON  HYDROLYSIS  RATE 

These  experiments  were  performed  at  a  constant  temperature  and  droplet  sizes  of  0.05,  0.2 
and  2.0  pL  with  40  pL  water  on  concrete  CB1-37C;  the  hydrolysis  rate  was  not  affected  by  the  drop 
size  (Figure  6a). 

ACTIVATION  ENERGY  FOR  DMMP  HYDROLYSIS  ON  CONCRETE  CHUNKS 

The  activation  energy  was  calculated  using  the  Arrhenius  equation.  These  two  equations  are 
generally  applied  to  homogeneous  reactions  in  solution.  To  the  extent  that  the  reactants  are  liquids 
within  the  pores  of  the  concrete,  these  equations  may  be  applied  here.  The  equations  used  were:4 
k  =  A  exp(-Ea/RT)  and  In k  =  InA  -  Ea/RT 

where  k  =  rate  constant,  A  =  preexponential  factor,  R=  gas  constant,  Ea  =  activation  energy,  and  T  = 
temperature  in  Kelvin. 

The  activation  energy,  Ea,  was  derived  from  the  slope  divided  by  the  gas  constant  R  from  a 
plot  of  Ink  vs  1/T„  and  the  preexponential  factor  A  was  calculated  from  the  intercept  of  the  line 
(Figure  6b,  Table  3). 


Figure  6.  a)  Half-life  for  DMMP  hydrolysis  at  40°C  as  a  function  of  droplet  size,  b)  Arrhenius  plot 
of  Ink  vs  1/T.  ♦=  CB1-27C,  24pL  water,  A=  CB1-37C,  8pL  water,  O  =  CB1-37C,  24pL  water, 

□  =  CB1-37C,  40pL  water. 


Table  1.  Summary  of  Activation  Energy  Parameters 


Concrete 

pL  water 

Ea  kJ  mol’1 

A,  M’V1 

-slope,  K 

CB1-37C 

8 

78.5 

3.1E10 

9446 

CB1-37C 

24 

79.1 

9.7E10 

9519 

CB1-37C 

40 

68.4 

2.2E9 

8230 

CB1-27C 

24 

77.6 

5.8E10 

9341 

DMMP  HYDROLYSIS  IN  AQUEOUS  SOLUTION 

The  purpose  of  these  experiments  was  to  simulate  the  hydrolysis  in  the  presence  of  heavy 
rainfall.  The  rate  generally  followed  that  of  DMMP  in  pH -buffered  solution,  and  was  base- 
catalyzed. 


Figure  7.  Plot  of  30°C  half-life  versus  pFI.  ♦=  Concrete,  40pL  water;  0=  Concrete,  8pL  water, 
■=  Aqueous  solution;  A=  Concrete  in  lmL  water. 

PRODUCT  ANALYSIS 

The  reaction  of  interest  was  the  hydrolysis  of  DMMP,  which  has  been  well  documented  in 
the  literature.  Methyl  methylphosphonic  acid,  MMP,  is  the  product;  the  second  hydrolysis  does 
not  occur.  The  reaction  has  been  shown  to  proceed  with  base  catalysis.5'6 

(CH30)2CH3P0  +  H20  -»  CH3OH  +  (CH30)CH3P00H 
DMMP  +  H20  CH3OH  +  MMP 

In  this  work  the  products  were  identified  using  31P,  ’H-coupled  and  ’H,  31P-coupled  NMR 
spectroscopy  (Figure  8). 
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Figure  8.  a)  31P,  ’H-coupled  NMR  spectrum  and  b)  *H,  31P-coupled  NMR  spectrum  of  DMMP  and 
MMP  in  aqueous  basic  solution. 


CONCLUSIONS 

Samples  of  concrete  from  different  sources  gave  different  rates  for  the  hydrolysis  of  DMMP. 
One  difference  found  between  the  concrete  from  various  sources  was  that  the  pH  imparted  to  water 
in  which  they  were  placed  differed;  it  ranged  from  8  to  1 1  pH  units. 

The  hydrolysis  of  DMMP  in  aqueous  solution  increased  with  increasing  pH;  its  hydrolysis  on 
concrete  chunks  followed  this  trend. 

For  concrete  CB1-37,  it  was  possible  to  quantify  the  hydrolysis  rate  in  terms  of  the 
temperature  and  amount  of  water  added. 

The  activation  energy  of  the  hydrolysis  of  DMMP  did  not  change  with  the  concrete  used,  but 
decreased  when  more  water  was  present. 


Changing  the  droplet  size  from  0.05  pL  to  2.0  pL  had  no  affect  on  the  DMMP  hydrolysis 

rate. 

This  work  will  be  extended  to  include  surface  area  and  pore  size  characterization  of  the 
concrete  samples,  and  more  types  of  concrete  and  aggregate  will  be  added.  Smaller  droplets  will  be 
used,  and  studies  of  VX  and  the  simulant  0,S-diethyl  phenylphosphonothioate  (DEPPT)  will  be 
performed. 


ACKNOWLEDGMENTS 

We  thank  Drs.  Vicky  Bevilacqua  and  Jeff  Rice  for  technical  assistance,  Wendel  Shuely  and 
Richard  O’Connor  for  the  concrete  samples,  Joe  Kiple,  Mark  Brickhouse  and  H.  Dupont  Durst  for 
supporting  the  program. 


REFERENCES 

1  Groenewold,  G.  S.;  Appelhans,  A.  D.;  Gresham,  G.  L.;  Olson,  J.  E.;  Jeffrey,  M.;  Weibel,  M.  J.  Am.  Soc. 
Mass  Spectrom.  2000, 11,  69-77 . 

2  Wagner,  G.  W;  O’Connor,  R.  J.;  Procell,  L.  R.  Langmuir  2001, 17(14),  4336-4341. 

3  Groenewold,  G.S.  et  al.  Envrion.  Sci.  Teclmol.  2002,  36,  4790-4794. 

4  Avery,  H.  E.  “Basic  Reaction  Kinetics  and  Mechanisms”,  pp  50-55,  1974,  The  Macmillan  Press  Ltd., 
London. 

5  Christol,  H.;  Levy,  M.;  Marty,  C.  J.  Organometal.  Chem.  1968,  12,  459-470. 

6  Christol,  H.;  Marty  C.  J.  Organometal.  Chem.  1968,  12,  471-478. 


